function p = psnr(x,y)
    %pour chaque pixel, moyenne de la différence au carré
    %10 * log10( erreur maximal (255² / calcul precedent)
    
    erreur = (double(x) - double(y)) .^ 2;
    [rows,cols] = size(x);
    moy_quadra = sum(erreur(:)) / (rows * cols);

    p = 10 * log10( 256^2 / moy_quadra);

end

